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Abstract 

The relative rates of l?-meson decays into J/ip and ip(2S) mesons are measured 
for the three decay modes in pp collisions recorded with the LHCb detector. The 
ratios of branching fractions (£>) are measured to be 

B B B [B+^f/^K+) = 0.594 ± 0.006 (stat) ±0.016 (syst) ±0.015 (i?^), 

T^i/ff? = 0.476 ± 0.014 (stat)± 0.010 (syst)± 0.012 (i?^), 

g fflffi/iy = 0.489 ± 0.026 (stat)± 0.021 (syst)± 0.012 (i?^,), 

where the third uncertainty is from the ratio of the ip(2S) and J/tp branching frac- 
tions to . 
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1 Introduction 



Decays of B mesons to two-body final states containing a charmonium resonance such 
as a J/ip or ip(2S) offer a powerful way of studying electroweak transitions. Such decays 
probe charmonium properties and play a role in the study of CP violation and mixing in 
the neutral B system (Tl. 

The relative branching fractions of B + , B° and B° s mesons into J/ip and ip(2S) mesons 
have previously been studied by both the CDF and DO collaborations [2}j4]. Since the 
current experimental results for the study of CP violation in B® mixing using the B® —> 
J/ijj (f> decay [5}|7] are statistically limited, it is important to establish other channels where 
this analysis can be done. One such channel is the B° s — > ip(2S)<p decay. 

In this paper, measurements of the ratios of the branching fractions of B mesons 
decaying to ip(2S)X and J/ipX are reported, where B denotes a B + , B° or B® meson 
(charge conjugate decays are implicitly included) and X denotes a K + , K*° or <p meson. 
The data were collected by the LHCb experiment in pp collisions at the centre-of-mass 
energy y/s = 7TeV during 2011 and correspond to an integrated luminosity of 0.37fb _1 . 



2 Detector description 

The LHCb detector M is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of b- and c-hadrons. The detector includes 
a high precision tracking system consisting of a silicon-strip vertex detector surround- 
ing the pp interaction region, a large-area silicon-strip detector located upstream of a 
dipole magnet with a bending power of about 4Tm, and three stations of silicon-strip 
detectors and straw drift-tubes placed downstream. The combined tracking system has 
a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and an impact parameter resolution of 20 |im for tracks with high transverse momentum. 
Data were taken with both magnet polarities to reduce systematic effects due to detector 
asymmetries. Charged hadrons are identified using two ring-imaging Cherenkov (RICH) 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, and electromagnetic and hadronic 
calorimeters. Muons are identified by a muon system composed of alternating layers of 
iron and multiwire proportional chambers. The trigger consists of a hardware stage based 
on information from the calorimeter and muon systems, followed by a software stage which 
applies a full event reconstruction. 

Events with a J/if> — > fi + fi~ final state are triggered using two hardware trig- 
ger decisions: the single-muon decision, which requires one muon candidate with a 
transverse momentum larger than 1.5 GeV/c, and the di-muon decision, which re- 
quires two muon candidates with transverse momenta p^ 1 and pt 2 satisfying the relation 
y/PT! ■ Pt 2 > 1-3 GeV/c. The di-muon trigger decision in the software trigger requires 
muon pairs of opposite charge with p^ > 500 MeV/c, forming a common vertex and with 
an invariant mass in excess of 2.9 GeV/c 2 . 
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3 Event selection 



In this analysis, the decays B + — y ipK + (B° — y ipK*°, B° s — y ip<f)) are reconstructed, where 
ij) represents ip(2S) or J/ip , reconstructed in the ip —y fi + fi~ decay modes. A K + (K*°, 0) 
candidate is added to the di-muon pair to form a B + (B°, B° s ) candidate. 

The starting point of the analysis is the reconstruction of either a J/ip or ip(2S) meson 
decaying into a di-muon pair. Candidates are formed from pairs of opposite sign tracks 
that both have a transverse momentum larger than 500MeV/c. Good reconstruction 
quality is assured by requiring the \ 2 P er degree of freedom of the track fit to satisfy 
X 2 /ndf < 5. Both tracks must be identified as muons. This is achieved by requiring the 
muon identification variable, the difference in logarithm of the likelihood of the muon and 
hadron hypotheses [9| provided by the muon detection system, to satisfy Alog£ M_?i > 
—5. The muons are required to form a common vertex of good quality (Xvtx < 20). 
The resulting di-muon candidate is required to have decay length significance from its 
associated primary vertex greater than 5 and have an invariant mass between 3020 and 
3135MeV/c 2 in the case of a J/tp candidate or between 3597 and 3730MeV/c 2 for a ip(2S) 
candidate. These correspond to [— 5a; 3a] windows around the nominal mass. The 
asymmetric window allows for the QED radiative tail. 

The selected J/ip and ip(2S) candidates are then combined with a K + , K*° or to 
create B meson candidates. Only the K*° — y K + tt~ and <fi —y K + K~ decay modes are 
considered. Pion-kaon separation is provided by the ring-imaging Cherenkov detectors. 
To identify kaons the difference in logarithm of the likelihood of the kaon and pion hy- 
potheses [9] is required to satisfy Alog£ x_7r > —5. In the case of pions the difference 
in logarithm of the likelihood of the pion and kaon hypotheses [9] is required to satisfy 
A log C K ^ K > —5. As in the case of muons, a cut is applied on the track % 2 /ndf provided 
by the track fit at 5. The kaons and pions are required to have a transverse momentum 
larger than 250 MeV/ c and to have an impact parameter significance with respect to any 
primary vertex larger than 2. In the B° channel, the mass of the kaon and pion system is 
required to be 842 < M^+n- < 942 MeV/ c 2 and in the B® channel the mass of the kaon 
pair is required to be 1010 < M K+K - < 1030MeV/c 2 . 

In addition, we require the decay time of the B candidate (cr) to be larger than 100 |xm 
to reduce the large combinatorial background from particles produced in the primary pp 
interaction. A global refit of the three-prong (four-prong) combination is performed with 
a primary vertex constraint and with the di-muon pair mass constrained to the nominal 



value [lO] using the Decay Tree Fit (DTF) procedure 11 . The reduced \ 2 °f this fit 
(Xprpp/ndf) is required to be less than 5, where the DTF algorithm takes into account 
the number of decay products to determine the number of degrees of freedom. The B + 
candidates, where a muon from the ip(2S) —y fi + n~ decay is reconstructed as both muon 
and kaon, are removed by requiring the angle between the same sign muon and kaon to 
be greater than 3 mrad. 
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Figure 1: Mass distributions of (a) B + -»• J/^K + , (b) B + -> ^(25)^+, (c) £° -)• J^K+w', 
(d) 5° -> ip(2S)K + n-, (e) 5° -»• J/^K+K- and (f) £° -»• V(2S)^ + i^~. The total fitted func- 
tion (solid) and the combinatorial background (dashed) are shown. The variation in resolution 
of the different modes is fully consistent with the energy released in the decays and in agreement 
with simulation. 



4 Measurement of N^^x/Nj/^x 

The mass distributions for selected candidates are shown in Fig. [T] The number of the 
B + — y ipK + candidates is estimated by performing an unbinned maximum likelihood fit. 
The same procedure is used to determine the number of the B° — > tpK + TT~ candidates in 
a 842 < Mjf+ r < 942 MeV/c 2 mass window and the number of the B® -y i[>K + K~ 
candidates in a 1010 < M^+k- < 1030MeV/c 2 mass window. The number of signal 



candidates is determined by fitting a double-sided Crystal Ball function 12,13 for signal 
together with an exponential function to model the background. The tail parameters of 
the Crystal Ball function are fixed to values determined from simulation. 

The B° mass distributions include the contributions from resonant decays (B° — y 
i/jK* ), non- resonant decays (B° — y ipK + 7T~) and combinatorial background. The contri- 
butions from resonant and non- resonant modes are separated with the sPlot technique [14] . 
The K + 7i~ invariant mass is used as a discriminating variable to unfold the B° mass dis- 
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Table 1: Summary of the signal yields for the six B modes considered and the ratios of the 
number of J/tp and ip{2S) decays: _/V total is the summed signal yield for resonant and non- 
resonant modes, 7\r non_res is the signal yield for non-resonant modes only and is the signal 
yield for resonant decays (through K*° or <p). The uncertainties are statistical only. 



B decay modes 


^ytotal 




jyres 


jvrres //V res 


B + J/ipK + 
B + ^(2S)K+ 


141, 769 ±410 
12, 154 ± 130 




141,769 ±410 
12, 154 ±130 


0.0857 ±0.0009 


B° J/tpK+n- 
B° ip(2S)K+ir- 


35, 770 ± 207 
2, 223 ±60 


1,253 ±30 
112 ±12 


34, 517 ±209 
2,111 ±61 


0.0612 ±0.0018 


B° s -> J/tfjK + K- 
B° s -> ^{2S)K + K~ 


7, 654 ± 92 
495 ± 25 


66 ± 13 
+1 


7, 588 ±93 
495 ± 25 


0.0652 ± 0.0034 



tribution of non-K*° K + tt~ combinations. A fit is then performed to the unfolded B° 
distribution, which contains both non-resonant B° — >■ ipK + ir~ decays and background, 
to determine the number of non-resonant decays. The final number of resonant decays is 
calculated by subtracting the number of non-resonant decays from the total number of 
decays. For the B® modes the number of non-resonant decays are obtained by a similar 
procedure using the K + K~ invariant mass as the discriminating variable. The signal 
yields and their ratios are summarized in Table [T} 



5 Efficiencies and systematic uncertainties 

The branching fraction ratio is calculated using 

B(B -» jj(2S)X) = N™ s)x ^ ej^ x B(J/i/j -> 
B(B -> J/^X) Nj; x X e H2S)x X Bty(2S) -> n+n~Y 1 J 

where N res is the number of signal candidates and e is the overall efficiency. 

The overall efficiency is the product of the geometrical acceptance of the detector, 
the combined reconstruction and selection efficiency, and the trigger efficiency. The effi- 
ciency ratio is estimated using simulation for all six modes. The simulation samples used 
are based on the Pythia 6.4 generator fTBI configured with the parameters detailed in 



Ref. 16 . Final state QED radiative corrections are included using the Photos pack- 



age 17 . The EvtGen 18 and Geant4 19 packages are used to generate hadron 
decays and simulate interactions in the detector, respectively. The digitized output is 
passed through a full simulation of both the hardware and software trigger and then 
reconstructed in the same way as the data. 

The overall efficiency ratio is 0.901 ±0.016, 1.011 ±0.014 and 0.994 ±0.014 for the B + , 
the B° and the B® channels respectively. Since the selection criteria for B — > J/ipX and 
B — > ip(2S)X decays are identical, the ratio of efficiencies is expected to be close to unity. 
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Table 2: Systematic uncertainties (in %) on the relative branching fractions. 



Source 


B + channel 


B° channel 


B° s channel 


non-resonant decays 


— 


1.5 


3.4 


data-simulation agreement 


1.7 


0.5 


2.0 


magnet polarity 


1.4 


0.6 


0.7 


finite simulation sample size 


0.3 


0.5 


0.6 


trigger 


1.1 


1.1 


1.1 


background shape 


0.6 


0.2 


0.2 


signal shape 


0.7 


0.8 


0.5 


angular distribution 




< 0.1 


0.6 


particle misidentification 


0.4 


< 0.1 


< 0.1 


Sum in quadrature 


2.7 


2.2 


4.3 



The deviation of the overall efficiency ratio from unity in the case of the B + — > tpK + 
decays is due to the difference between the spectra of muons for the J/ip and ip(2S) 
decays. For the B° and B s channels this difference is small. It has been checked that the 
behaviour of the efficiencies of all selection criteria is consistent in the data and simulation. 

Since the decay products in each of the pairs of channels considered have similar 
kinematics, most uncertainties cancel in the ratio. The different contributions to the sys- 
tematic uncertainties affecting this analysis are discussed in the following and summarized 
in Table H 

The dominant source of systematic uncertainty arises from the subtraction of the non- 
resonant components in the B° and the B° s decays. The non-resonant background is 
studied with two alternative methods. First, determining the number of B9^ — > ipK* Q {(j)) 
decays directly using the sPlot technique by unfolding and fitting the mass distri- 
bution of candidates containing genuine K*°((fi) resonances. Second, using the B^ mass 
distribution as the discriminating variable to unfold the K + ir~(K + K~) mass distribution 
of genuine B^ candidates and fitting this distribution to determine the number of non- 
resonant decays. The corresponding uncertainties are found to be 1.5% in the B° channel 
and 3.4% in the B° s channel. 

The other important source of uncertainty arises from the estimation of the efficiencies 
due to the potential disagreement between data and simulation. This is studied by varying 
the selection criteria in data and simulation. The corresponding uncertainties are found 
to be 1.7% in the B + channel, 0.5% in the B° channel and 2.0% in the B° s channel. The 
observed difference in the efficiency ratios for the two magnet polarities is conservatively 
taken as an estimate of the systematic uncertainty. This is 1.4% in the B + channel, 0.6% 
in the B° channel and 0.7% in the B° s channel. 

The trigger is highly efficient in selecting B meson decays with two muons in the final 
state. For this analysis the di-muon pair is required to trigger the event. Differences in 
the trigger efficiency between data and simulation are studied in the data using events 
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which were triggered independently on the di-muon pair 20 . Based on these studies, an 
uncertainty of 1.1% is assigned. 

A further uncertainty arises from the imperfect knowledge of the shape of the signal 
and background in the B meson mass distribution. To estimate this effect, a linear 
and a quadratic function are considered as alternative models for the background mass 
distribution. In addition, a double Gaussian shape and a sum of double-sided Crystal Ball 
and Gaussian shapes are used as alternative models for the signal shape. The maximum 
observed change in the ratio of yields in the ip(2S) and J/ip modes is taken as systematic 
uncertainty. 

The central value of the relative efficiency is determined by assuming that the angular 
distribution of the B — > %p(2S)X decay is the same as that of the B — > J/ipX. The 
systematic uncertainty due to the unknown polarization of the i/j(2S) in the B meson 
decays is estimated as follows. The simulation samples were re-weighted to match the 
angular distributions found from the data and the relative efficiency was recalculated. 
The difference between the baseline analysis and the re-weighted simulation is taken as 
the systematic uncertainty, as shown in Table [2j 

Finally, the uncertainty due to potential contribution from the Cabibbo-suppressed 
mode with a 7r misidentified as K is found to be 0.4% in the B + channel and negligible 
in the B° and B® channels. The uncertainty due to the cross-feed between B° and B® 
channels with a 7r misidentified as K (or a K misidentified as it) is negligible. 



6 Results 

Since the di-electron branching fractions are measured more precisely than those of 
the di-muon decay modes, we assume lepton universality and take = B(J/ip — > 
fi + Li-)/B(4>(2S) n+yr) = B{J/ip -» e + e-)/B(if)(2S) e + e~) = 7.69 ± 0.19 (To] . 
The results are combined using Eq. [T] to give 

B ^B?*ffrK+) = 11594 ± 0.006 (stat) ±0.016 (syst) ±0.015 (i^), 

%*&*W%P> ) = 0.476 ±0.014 (stat)± 0.010 (syst)± 0.012 (^), 

g gflg^!wff = °- 489 ± 0.026 (stat)± 0.021 (syst)± 0.012 (^), 

where the first uncertainty is statistical, the second is systematic and the third is the 
uncertainty on the value 10 . 



The resulting branching fraction ratios are compatible with, but significantly more 
precise than, the current world averages of B(B + — > ip{2S)K + ) /B(B + — > J/ipK + ) = 
0.60 ± 0.07 and B(B° S -> ip(2S)(f>) /B(B° S -> J/ij}(j>) = 0.53 ± 0.10 |10| and the CDF result 
of B(B° -> tP(2S)K*°)/B(B° ->■ J/ipK* ) = 0.515 ± 0.113 ± 0.052 |2|. The B° s ->■ ip(2S)<p 
decay is particulary interesting since, with more data, it can be used for the measurement 
of CP violation in B° s mixing. 
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